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Abstract. We discuss the occurrence and properties of exact singular anyonic wave functions describing 
stationary states of two identical charged particles moving on a plane and under the influence of a perpen- 
dicular uniform magnetic field. 

PACS. 03.65.Ge Solutions of wave equations: bound states - 05.30.Pr Fractional statistics systems - 
71.10.Pm Fermions in reduced dimensions - 73.43.-f Quantum Hall effects. 



1 Introduction 

The general principles of quantum mechanics requires that 
localized stationary states of a system be represented by 
square integrable functions so as to be consistent with the 
probabilistic interpretation of Born. This is in particular 
the situation for the motion in a potential well rising in- 
definitely at infinity and at a singular point localized, say 
at the origin. Here we are concerned with this condition 
in two space dimensions. Let us suppose that the wave 
function \P behaves as r s u(r) near the origin where u(r) is 
a well-behaved function of r, the distance from the origin, 
for r — > 0. Then for J \tif\ 2 d?r, to be convergent in a neigh- 
borhood of the origin, one should have s > — 1. There is 
however a second condition stating that the origin as a 
singular point should not be a source nor sink for matter 
The outward flow of the probability flux J through a 
small sphere of radius r must vanish as r — ► 0: 
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m being the particle mass. This puts a stronger bound on 
s: i.e. s > — h. However if u is a real function this condition 
is always satisfied as it is in the case we consider now. 

In most soluble two-dimensional problems presenting 
cylindrical symmetric quantization of angular momentum 
to values hi leads to the condition s 2 = I 2 or s = ±\l\ 
equivalently. For the wave function to be single- valued one 
imposes I 6 Z and the choice s = is excluded since 
it is inconsistent with s > —1. But with the discovery 
of anyonic behavior for wave functions of identical planar 
particles by Leinaas & Myrheim I as quantized rela- 
tive angular momentum of two particles is a priori a real 
number: I E R. Then it is possible to consider the sec- 
ond option s = so long as s > — 1. This possibility 
which has been pointed out a decade ago by Grundberg et 



al H in a general context leads to singular anyonic wave 
functions. 

In this short note, we present examples of such func- 
tions occurring in the relative motion of two identical pla- 
nar particles of mass m and charge e in a perpendicular 
uniform magnetic field B. First we show how to construct 
the wave functions in terms of Biconfluent Heun Polyno- 
mials and how quantization conditions are imposed. In 
the next section we discuss the structure of the low-lying 
singular wave functions. Finally we close the paper with 
some comments and physical interpretations on this un- 
usual behavior. 



2 Construction of the wave functions 

As the Coulomb repulsion between the particles has ro- 
tational symmetry, we use cylindrical coordinates and the 
symmetric gauge to describe the magnetic field. The Hamil- 
tonian of this relative motion reads: 
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Here, r is the particle separation, 8 its inclination with 
respect to a reference axis, p r and pe are the conjugate 
momenta and K — -j— in S.I. units. 

The wave functions of stationary states ^(r, 6) are so- 
lutions of the Schrodinger equation HlP = E\P with energy 
E. As pg is conserved and in view of possible exotic statis- 
tics in two dimensions we may set &{r,9) ~ e ue Ri(r) with 
I £ R. To put forward the competition between magnetic 
length Ib and the Bohr radius Iq defined by: 
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we recast the radial Schrodinger equation for Ri (r) in its 
dimcnsionlcss form: 
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using the scaled variables £, e and C given by: 
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/From cq. (0) , it can be seen that the asymptotic form 

of Ri(£) is e~2~. We set Ri(£) — Vi(O e ^~ w i tn the as- 
sumption that vi(£) is of milder growth, at £ — > oo, than 

e 2 . «/(£) fulfills now an "associated" Biconfluent Hcun 
equation: 

+ (1 - 2£ 2 )» J ' - {C + 2(e - / - - -t;, = (5) 

with Maroni's canonical parameters a = /3 = 0, 
7 = 2(e — Z) and <5 = 2C Q. We seek now a power series 
solution of the form = EriLo a ™£ s+ "- Substitution 

of this series into eq.(y) yields the recursion relations: 
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Assuming «o ^ 0, we have s 2 = I 2 or s = <j[ 
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we recover the 



a = ± . Now posing a„ = (1+2g)nra! 
canonical recursion of the Biconfluent Hcun Function 
N(2s,0,2(e-l),2C;0 |: 

-2(n + l)(n + 1 + 2s)(e - n - 1 - I - s)A^ , (7) 

with = 1 and A^" 1 = C; the parameter s may be +\l\ 
or — \l\ and label regular or singular wave functions near 
the origin. 

As AT(2s, 0, 2(e - l),2C; f) grows as {§ it is neces- 
sary to cut its defining power series down to a polynomial 
by imposing two conditions on the recursion relation (Q): 

— energy quantization e = e n = (n + 1 + I + s), 

- and A ( °l 1 (C,l)=0. 

In general for given B (or C) the angular momentum 
I = Ik is the fc" 1 solution of an algebraic equation of order 
n for n odd and of order (n — 1) for rt even. This has 
been already pointed out in || . The main consequence is 
that the magnetic field fixes I and thus gives rise to the 
expected anyonic behavior. Note that the effective radial 
potential becomes also modified for each eigenstate since 
it depends on I. 



Globally the wave function of a stationary state reads: 

(8) 



where n n (£, C) is a special Biconfluent Hcun polynomial, 
obtained by truncating the power series. 

As well-behaved wave functions with a = + have been 
already discussed in Q, we shall be concerned here with 

the low- lying with \h\ < 1< The corresponding 

energy eigenvalues are e^j = n + 1 + Ik — \h\- In the 
following we shall assume that C ^ (no infinitely strong 
magnetic field). 



3 Properties of the wave functions 

They appear to be similar for pairs of the quantum num- 
ber n. This is essentially due to the fact that the conditions 
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for n even lead to a pair of 



polynomials of same order in C 2 since A^}^ contains an 
overall factor C. For a = — , we shall study the variation 
of C 2 as a function of I in the interval — 1 < I < 1. Two 
possible behaviors of the radial probability density for the 
relative separation r arise: divergent for ^ < |Z| < 1 and 
convergent for < \l\ < | in the neighborhood of zero. 
This is the only physical relevant quantity since rotational 
symmetry yields a uniform probability density for the in- 
clination angle 6. 

The first group of states is related to n = 1,2. The 
vanishing of A^h = yields linear relations between 
C 2 and I, For n = 1,2, the equations are respectively 
C 2 = 2(1 - 2|/|) and C 2 = 4(3 - 4|Z|). This is shown in 
Fig.@. We have given the probability distribution of states 
r\ip~ j (9, r)\ 2 in Fig.|| and Fig.|| as function of r. 
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To save space, normalization factors with involved expres- 
sions are not given but are implicitly taken into account. 
In particular, we see that at some value of C < 2 for 
n = 2 the radial probability distribution starts to develop 
an integrable singularity at the origin. 

The next group of states is given by n = 3,4. The 
conditions fixing I in terms of C are A{ \C,1) = and 
A\'{C, I) — 0, which correspond to the equations C 2 = 
20(1 - |Z|) ± ^/64|Z| 2 -28|Z| + 73 and C 2 = (50 - 40|/|) ± 

6^16|Z| 2 -40|Z| + 33. Since C = is discarded we obtain 
quadratic equations which yield two values If., k = 1,2. In 
Fig.||,(rcsp. p|) we plot C 2 in terms of I for n — 3, (respA): 
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n=l and n=2 with s=— 1 
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Fig. 1. C 2 versus I for n = 1 (broken line) and n — 2 (contin- 
uous line) 



1 value of C 
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Fig. 2. Probability distribution P{ (r, C = 1) for n = 1 



2 values of C 
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Fig. 3. Probability distributions P 2 ( r i C) f° r n = 2, with 
C = 2^2 (cont inuous line) and C = v2 (broken line). 
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Fig. 4. C versus of I respectively for n = 3, the upper (resp. 
lower) branch correspond to + (resp. — ) in expression of C 2 . 

We observe that for certain values of the magnetic held 
such that C 2 is in a gap: 

- (20 - 2^73) < C 2 < 6 for n = 3 

- (50 - 6\/33) < C 2 < 28 for n = 4 

there can be no states 3^ j (0, r) for n = 3, 4. The behavior 
of the probability distributions depends on whether |Z| is 
larger or smaller than ~. 

For n = 3we may have 3 types of curves (see Fig.||): 
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Fig. 5. Probability distributions P 3 (r, C) for n — 3 
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Fig. 6. C 2 versus of I respectively for n — 4, the upper (resp. 
lower) branch correspond to + (resp. — ) in expression of C 2 

- Pt\r, C = 5) (broken line) and p!f ] (r, (7=1) (con- 
tinuous line) have smooth behavior at the origin. 

- Pt\r,C = 2V2) (dotted curve, with magnified part 
in lower inset) is divergent but integrable at the origin. 

For n = 4, we present 2 pairs of curves in Fig.^: 

- P ( f\r,C = 9) (continuous line) and P^\r,C = 3) 
(thin broken line) which go smoothly to zero at the 
origin 

- Pt\r,C = 6) (bold broken line) and Pt\r,C = 2) 
(dotted line) which diverge at r = but remain inte- 
grable in this neighborhood. This curves are magnified 
in insets in the lower part of Fig.|^. 

The P 3 4 (r, C) are computed with the expression of 
the wave function in eq.(^). Since the form of the spe- 
cial Biconfluent Heun Polynomials .773,4 (£, C) are quite in- 
volved, we shall not give here the explicit form of P3 ^ (r, C) 

Higher stationary states come also in pairs n — 5, 6 and 
n = 7, 8 with algebraic relations of third and fourth order 
in C 2 (and/or I). The angular momentum I will have re- 
spectively three and four values and a more complicated 
structure in terms of magnetic field arises. General fea- 



4 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 



4 values of C 
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Fig. 7. Probability distributions P\ (r, C) for n = 4 



tures concerning two classes of smooth and singular be- 
havior at the origin however remain. We shall discuss them 
in a more extensive study later on. Beyond these values of 
n, one can only tackle the problem numerically since the 
algebraic equations of higher order do not posses solutions 
in close form. 



n 


B 


C 


h 




n 
la 


1 


940215 


1 





3 


0.3409025 


2 


470108 


V2 


0.207107 


1.75 


0.5 




117527 


2V2 





0.625 


0.380453 


3 


940215 


1 





6.60434 


0.215634 




37608.6 


5 





0.263969 


3.84573 




117527 


2V2 


0.322087 


0.545623 


1.99029 


4 


235054 


2 


0.1032474 


1.9443 


0.666414 




104468 


3 





0.965468 


1.07414 




26117.1 


6 


0.386699 


0.184408 


5.4233 




11607.6 


9 





0.120639 


9.1545 



4 Comments and conclusion 

The presence of fractional power singularities of the wave 
functions at the origin of coordinates calls for some com- 
ments. The rise of the probability of finding a separation 
r between two identical particles within dr suggests a ten- 
dency for the particles to cluster together. Since 
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and as < A < 1, 



integration of r from to any value r = a close to the 
origin is finite. This means that the probability of find- 
ing any r in the interval [0, a] is always finite, despite the 
apparent divergence. Moreover the clustering tendency in 
the presence of classical Coulomb repulsion is somewhat 
counter-intuitive. To reconcile this behavior with quan- 
tum mechanics we have given in Table. [I] the values of the 
zeros ro of the probability distributions nearest to zero, 
from which start divergent behaviors for various n and Ik 
as well as the locations r% below which classical motion 
is not possible in corresponding effective potentials. For 
the four low lying n, we see that always rg < r\. Hence 
the divergence occurs only in forbidden classical regions 
of the motion, thus not directly observable. Anyhow, such 
situations occur also at tremendously high B-ficld, of the 
order of 10 5 Tesla. 

Finally, it is amusing to note that there is an analog 
of this situation in electrostatics. The electric field, near 
the edge of a two-dimensional wedge of opening angle 



is the analog of ip, 
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") and behaves as 



be- 



ing the distance from the edge |qj. Thus it diverges for 
7r < 6 < 2ir. Now the electric field energy density is the 
analog of P~ (r, C) , but here it has always a nice behavior 
r 2 !" 1 for < 6 < 2ir. 
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Table 1. Table of closest zeros ro of low- lying probability 
distributions with given C, energy levels and the locations of 
the classically limiting points n at these values of the energy 



(-) 



n + 1 — 2\l\). Values of the magnetic B in Tesla, corre- 



sponding to chosen values of C are listed in the second column. 



